Solution of the asymmetric double sine-Gordon equation 
Nan- Hong Kuo 1 , Sujit Sarkar 2 , C. D. Hu 1 ' 3 
department of physics, National Taiwan University, Taipei, Taiwan, R.O.C. 
2 Poornaprajna Institute of Scientific Research, 4 Sadashivanagar, 
Bangalore-560 080, India. 
3 Center for Theoretical Sciences, National Taiwan University, Taipei, Taiwan, 

R.O.C. 
Abstract 

We present solutions of asymmetric double sine-Gordon equation (DSGE) 
of an infinite system based on Mobius transformation and numerical exercise. 
This method is able to give the forms of the solutions for all the region on the 
ip — r\ parameter plane where ip is an additional phase and r\ is the ratio of 
the magnitudes of two sine terms. We are able to show how the deconfinement 
occurs near ip — (1/2 + n)7r and if = nit. and also find the solution for all 
values of <p. We predict different kind of solutions and transitions among them 
in different parts of the parameter space of this equation. 
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0.0.1 1. Introduction: 



The sine-Gordon equation attracted much interest of physicists [1,2]. In quan- 
tum field theory it is a rare example of integrable system which can be a starting 
point of developing non-perturbative theory. It also has plenty of application in 
condensed matter systems[3,4] and nonlinear optics[5]. However, the Lagrangian 
of a realistic physical system often gives a more complicated equation of motion 
than the sine-Gordon equation. For example, a quantum spin chain is mapped 
into a Lagrangian with several potential terms [6]. Systems with nonlinear op- 
tical properties also give rise to more complicated wave equations [7]. Thus a 
more complete model is desirable. This leads to the double sine-Gordon equation 
(DSGE) and loses the integrability of the sine-Gordon equation. But this "flaw" 
also provides a test ground for various analysis [8- 11] and perturbative or non- 
perturbative methods such as form factor perturbation theory[12, 13], .truncated 
conformal space approach [13, 14] and semiclassical approach[15]. For a standard 
DSGE, exact solutions can be found[8,10]. But for a asymmetric DSGE 

Ott- 6 xx + sin(8 + tp) +2r?sin 20 = 0. (1) 

with ip being an additional phase, exact solution has eluded the effort of re- 
searchers. It is interesting not only because the state equations of a strongly 
correlated electron system [16], a quantum spin chain with external field or that 
of a spin pump[17,18] can be mapped into eq. (1) but from a purely theoretical 
point of view, it contain rich physics as we shall see. 

The main difficulty of DSGE is that it is non-intcgrable. A rigorous ana- 
lytical form of the solutions of asymmetric DSGE is important because it can 
provide deeper insight to understand the physical systems. It is also well known 
that the perturbation theory is not applicable in one-dimension because there 
is no quasi-particle excitation alike in higher dimension. Various approximation 
methods, perturbative or non-perturbative, sometimes give results leading to 
different physical pictures[14]. For this reasons, we would not like to solve this 
equation with above methods. We propose an insightful analytical form of the 
solution which can serve as an anchor to numerical analysis or approximations. 
It can also be a springboard to study quantum fluctuations. 

The plan of the manuscript is the following: In Section 2, we present previ- 
ously existed classical solutions of DSGE depend on the potential coefficient i]. 
We also present explicitly shortcomings of these analytical solutions. Section 3 
contains the major part of our work. We present mathematical solution we find 
with Mobius transformation and see how parameters change with 77 and <p. In 
Section 4, We present results and discussions. 

0.0.2 2. Classical solutions of Double Sine-Gordon Equation 

In this section, we briefly review the previously existed solutions in the literature, 
potential and energy of the double sine-Gordon equation (DSGE) with <p = 
in order to establish the notations and set the stage for developing our method. 
We shall use the notations of ref. 4 and also point out the shortcoming of that 
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work. As mentioned in ref. 4, the Hamiltonian of the DSGE can be viewed as 
a chain of physical pendulums joined by torsion springs 

/2 p 
dx{±£ + -6l-V [cos6 + r]cos26}} (2) 
Zl z 

where 6 is the angular coordinate, p§ is the conjugate momentum, I is the 
moment of inertia, T is the torsional constant, and Vb is the external potential. 
The DSGE which we are interested in (eq. (1).) can be obtained with rescalings 
of the time-space coordinates. In order to calculate the total energy, one can 
Integrate both sides of eq. (1) and get 

1 df) 

-(^) 2 +cos(8 + <p)+r,cos(2e) = S (3) 

2 as 

where s = 7(2; — vt) , 7 = l/Vl — v 2 and v is the velocity of soliton. The 
negative sum of the last two terms can be viewed as the potential 

V(9;r),ip) = -cos(6 + ip) - r/cos(26») (4) 

and so S is called the " action" . We can show that 



S = -V min . (5) 

The solutions of kinks and bubbles (see ref. 4 and/or eqs. (28), (35) and (37) 
below) approach constant when s — > ±00. Hence, d6 / ds\ s= ± oc = and S = 
-V(0(±oo)). Since {d6/ds) 2 /2 = V(9) - V(9(±oo)) > for any s, we conclude 
that S = —V(0(±oo)) = — Vmin. More specifically, whether the minimum is an 
absolute minimum or a relative minimum depends on what type of the solution 
is. 

Now we would like to calculate energy for all ip. The calculation of energy 
can be performed with the a method similar to that of ref. 2. We get from eq. 
3. that 

d9 

— = ^2[S - cos(6» + ip)- 77 cos 26l]. (6) 
ds 

We define 

/°° f) 2 9 2 f°° 9 2 (1 -4- 7j 2 W f 2n 

fe + %]dx = (1+V 2 h / -fds = { + n / y/S-COs(e + <p)-TICOB20<W 
-00 2 Z J-00 2 v2 Jo 

(7) 

and 

/OO -j^ poo 

[-cos(6 + <p)- r] cos 28}dx =- / [-cos(6 + <p) - 77 cos 26>]ds. (8) 
-00 T J — 00 

so that H = V s + V p with some scaling. From eq. (3) we found that 

q . fi |s=i 1 roc n2 
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But the first term on the right hand side of eq. (8) diverges if we consider 
an infinite system (L — > oo). Thus we subtract this trivial infinity from the 
Hamiltonian in eq. (2). This is equivalent to shift the system to a new energy 
zero. In view of eq. (2), we get 

/oo rp. 
-ids. (10) 

While solving for solutions, we can obtain numerical values of S. Hence, the 
energy at every tp can be calculated. 

An analysis of the potential term is in order. It helps us to understand how 
state evolves with respect to the phase tp. The potential of double sine-Gordon 
equation has the symmetry 

V(0;ri,<p=^) = V(6;-ri,<p = 0) (11) 

where 9 = 9 + n/2. This implies that the solution of tp = n/2 has the same 
form as that of ip = and —77. However, interestingly, the solutions can have 
completely different forms for ip = or ip = n/2 when |7j| > \ (see below). 
In view of eq. (4), we have 

dV 

— =sm(9 + p) + 2r 1 sm(29) 
(19 

= sin 9 (cos ip + cot 9 sin ip + 4t7 cos 9) 
= cos 0(sin ip + tan 9 cos tp + Arj sin 9) . 

The point 9 = cos _1 (— 1/47/) is one of the absolute maxima only if pi = 
and 77 < —1/4, and 9 = sin _1 (— 1/4^) is one of the absolute minima only if 
<p = tt/2 and 77 > 1/4. Moreover, for 77 < —1/4 the relative minima are at 
9 = n/2 + 2nir and for 77 > 1/4 the relative maxima are at 6 = 2nn. The shape 
of the potential has critical influence on the form of the solution. We show 
their features in Figures la~lc versus 9 for rj = —0.35, 77 = 0.15 and 77 = 0.35 
respectively. Notably, for (77 = —0.35, tp = n/2) and (77 = 0.35, p = 0) there are 
relative minima and for (77 = —0.35, ip = 0) and (77 = 0.35, p = n/2) there are 
relative maxima. Nevertheless, for 77 = 0.15, there is no relative extrcmum. In 
general, there is a region where no relative extrcmum exists. In the regions where 
relative extremum exist, a so called bubble solution can be found. Another 
symmetry can be seen by making transformation & = 9 + <p, eq. (4) becomes 
V(6;r], tp) = — cos 9' — rjcos(29' — 2p). This shows that there is a period of n 
for the variation of ip. 

Specifically, when tp = and p — n/2 the analytical forms of solution and 
energy can be obtained. In order to set up our analysis, in the following we give 
a summary of the solutions and energies in the case tp = 0. 

Case 1 77 < — i 

In this case the absolute minima are at 9 a b s . m in = Pa + 2n7r where 

po = arccos( — ). (12) 
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From above, we have S = —V(6 ln i n ) = —r/ — 1 / 877. The absolute maximum 
of V{9) are located at # b s . max = (2n + l)ir while the relative maximum are 
located at 9 reLmax = 2mr with V(9 abs . max ) = 1 - 77 and V(0 rei . max ) = -1 - 77. 
There are two kinds of traveling kinks: 



9> = 2arctan[±J^coth J±^ s) ], (13) 



9< = 2 arctan[±^ ^ tanh(^ (14) 

where the superscripts > and < denote the large kink and small kink respec- 
tively. We discuss their energies separately. 

(A) Large kink: We found from eq. (13) that e [tp , 2n — ip ], and it must 
vary cross one of the absolute maxima. In view of eqs. (7) and (10), 

-^s = \/2 7 / J-t} cos6 -r)cos26d6 = \/2 7 / J-2r](cos9 + — ) 2 

-00 2 J Vo V 877 J va V 477 

:(Vl6r? 2 - l + 7T-^o) (15) 



(B) Small kink: From cq. (14) we found that 9 G {—Lpo,<fo], an d it must 
contain one of the relative maxima 



(H) = 2 7 f 6 -^ds = V27 / ^ J -rj - — - cosO - r ) cos26d6 = -^={^/l&7f - l-(p ) 
J-00 2 J V0 V 877 

(16) 

Here, we have corrected the errors in eqs. (3. 10) and (3.11) of ref. 4. 
Case 2 |t?| < \ 

There is only one type of basic kink solution in this case: 



9> = 2 arctan[± y/l + Ar/csch^l + 4r?s)] (17) 

The minimum of V{9) are located at 9 m i n — 2mr with V(0 m in) = — 1 — 77 and 
the maximum of V{9) are located at max = (2n + l)ir with V(9 max ) = 1 — 77. 
There is no relative maximum or minimum. We have S — — V(9 m [ n ) = 1 + 77 
and from eqs. (6) and (9), 

/oo n2 r-2ir 
-fds = \/2 7 / V I + 77 -cos 9 -77 cos29d9 (18) 
-00 2 Jo 

The integral also depends on wether 77 is larger or smaller than 0. 

(A) i > n > 



J -co A y/V 



5 



(B) > T) > 
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2-f arcsin(2- v /^) 



(H) = 2 7 / -fds = 4 7v ^+l + -i ^ v " (20) 

Here we have corrected an error in eq. (3.7) of ref. 4. 
Case 3 t] > j 

There are two kinds of traveling kinks: 



6 > = 2 arctanti-y/l + At]csch{y/l + 4r?s)], (21) 

and 



6 B = 2 a.rrta.n[+ — COSh(-v/4?7 - Is)] (22) 

V4?7 - 1 

where the superscript B denotes the bubble solution. The absolute minimum of 
V(9) are located at 9 a b s . min = 2ti7t with V{6 a bs. min) = — 1 — V and the maximum 
oiV(9) are located at 9 maK = arccos(— l/4r;) + 2n7r with V(9 a f )s . niax ) = 1/877 + 77. 
The energy of the large kink is the same as that in case 2: 



(B) - 2 T £ f d s - 4 TV ^TT + j^gj (23) 

The other kind of solution is the bubble solution. It extends from one relative 
minimum to another. These minima are at 6V e ; m m — > (2n + l)7r as s — > ±00. In 
this case, 5" = -V(0 re j,min) = 77 - 1 

/oo n2 i-tv 
-fds = 2V27 / y/r]-l-cos6-r]cos26d6 
-00 2 J2arctan(l/\/4»7-l) 

(24) 



= 47^/477 - 1 - 



4 7 ln(2^/77 + V4^7~~T) 



Note the upper bound and lower bound of the integral. We divide the bubble 
into two equal halves. The upper bound tt is the relative minimum while the 
lower bound is the middle point of the bubble. It comes from 9 B (s = 0) = 
2arctan[±(477 - 1)~ 1/2 cosh(^/4?7 - ls)]| s = = ±2 arctan( l/y/$rj - 1). Here we 
also corrected the errors in eq.(3.7) and (3.9) of ref. 4. 

Further insight can be gained by applying eq. (11). For example, if we start 
from tp = and 77 > j, the solutions are the large kink in eq. (21) and bubble 
in eq. (22). When (p is changed adiabatically into ir/2, the solutions become 
the large and small kinks in eqs. (13) and (14): 



9> = 2arctan[± W-i^ coth( 



9< = 2arctan[±4/4n tanh ( 



<Mv\ 


- 1 


4|r?| 


i 

+ 1 


U\v\ 


- 1 
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with the additional term —tt/2 coming from the difference between 9 and 9. 
Therefore, the solutions can have quite different forms as ip varies. The interest- 
ing question is whether the solution evolve smoothly or they change abruptly. 

0.0.3 3. Solutions in general 

In previous studies of asymmetric DSGE, the solution can be found when ip = 
or ip = ir/2. Here, we present a method which enables us to find solutions for 
any value of ip. We propose that the solution of cq. (1) in general has the form 

= 2arctan[/(s)]. (25) 

Then after substitution, we have the following equation 

2(f) 2 = (S+cos(<p)-r 1 )f 4 +2sm(v)f+(2S+6r 1 )f+2sm(<p)f+(E-co S (v)-r 1 ) 
as 

(26) 

Above equation is similar to the differential equation of Jacobi elliptic functions 
(JEF) or hyperbolic functions except for the terms with the odd power of f(s). 
However, for JEF, there is the "Mobius transformation" to change eq. (26) 
into the standard form. The details is given in Appendix. By letting f(s) = 
(ag(s) + b)/(cg(s) + d) and choosing suitable coefficients, a, b, c, and d, it is 
possible to obtain the following form from eq. (26) 

(J) 2 = a i9 4 + a 2 g 2 + a . (27) 

In an infinite system with arbitrary value of ip, it is reasonable to use for g(s) 
the hyperbolic functions which JEFs approach in the limit of modulus k — ► 1: 

= asinhM + 6 
J y 1 csinh(rs) + d y ' 

where r is a constant. One may also use the other hyperbolic function for this 
construction. Substituting eqs. (25) and (28) into cq. (26) and requiring the 
same scaling 

(ad - be) 2 = 1 (29) 
we have found the following equations by comparing the powers of sinh function 

(a 2 + c 2 ) 2 5'+(a 4 -c 4 )cos( / 9 + 2ac(a 2 + c 2 )sin^+(-a 4 + 6a 2 c 2 -c 4 )?7 = 0, (30a) 



A(ab + cd)(a 2 + c 2 )S + A(a 3 b - c 3 d) cos p 

+ 2(a 3 d + 3a 2 bc + 3ac 2 d + be 3 ) sin <p 

+ [-Aa 3 b + 6(2a 2 cd + 2a6c 2 ) - Ac 3 d]r/ = 0, (30b) 
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[6a 2 b 2 + 2(a 2 d 2 + Aabcd + b 2 c 2 ) + 6c 2 d 2 ]S 

+ (6a 2 b 2 - 6c 2 d 2 )cos<p> 
+ 6(a 2 bd + ab 2 c + acd 2 + bc 2 d) sin ip 
+ 6(-a 2 b 2 + a 2 d 2 + Aabcd + b 2 c 2 - c 2 d 2 )r] = 2r 2 (30c) 

A{ab + cd){b 2 + d 2 )S + A{b 3 a - d 3 c) cos p+ 

2(b 3 c + 3b 2 ad + 3bd 2 c + ad 3 ) sin tp+ 

[-4b 3 a + 6{2d 2 ab + 2cdb 2 ) - M 3 c]r] = (30d) 

(b 2 +d 2 ) 2 S+(b 4 -d i )cos(p + 2bd(b 2 +d 2 )sm(p+(-b 4 + 6b 2 d 2 -d 4 )r ] = 2r 2 (30c) 

The simultaneous algebraic equations are solved to give the coefficients a, 6, c, 
d, r and S. Thus, we can find solutions for any values of i] and <p except for 
certain special cases with the form in eq. (28). In Table 1, 2 and 3 we give the 
values of the parameters for several values of r\ and (p. The crosses in Table 3 
indicate that eq. (28) cannot give any solution under these conditions. 

The special case for which the form in eq. (28) cannot give any solution is 
well-defined and very interesting, it occurs at = tt/2 and rj > 1/4. We also 
explain the cause of this forbidden solution by starting with the symmetry of 
the DSGE with additional phase. If one substitutes it — ip for ip, then he can 
replace with tt — 9 and eq. (1) retains its original form. In view of eqs. (25) 
and (28), we have the following symmetry: 

ip — > 7T — p 

a — > — c 
b -» d 

c -» -a . (31) 
d -» b 
r — > r 
5 -» 5 

In other words, when <p is changed into tt — ip, the solution becomes tt — 9. 
Suppose we have, and indeed we have found, real solution of eqs. (30) for 
T] > 1/4 with a, b, c and d being continuous with varying ip, then we must have 
a = — c and b = d at ip = n/2. Substituting these into eqs. (30), we get 

a 4 (S - 1 + r,) = 0, (32a) 

b\S+l + ^) = r ^, (32c) 
along with the scaling restriction from eq. (29) 

Aa 2 b 2 = 1. (33) 
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Eq. (33) means a^O, and so cq. (32a) implies S + r] = 1. Inserting these into 
cqs. (32b) and (32c), we get 



r 2 



1. (34) 



I-477 

Incompatibility results from eq. (32) if r\ > 1/4. 

On the other hand, we can find a solution of distinctively different form at 
r) > 1/4 and tp = ir/2. Instead of g(s) = sinh(rs), we now have g(s) = tanh(rs) 
or g(s) — coth(rs) 

, . _ asinh(rs) + bcosh(ra) 
csinh(rs) + rfcosh(rs) 

The algebraic equations are 
(a 2 + c 2 ) 2 S+ (a 4 -c 4 ) cos p + 2ac(a 2 + c 2 ) sin p + {-a 4 + Qa 2 c 2 -c i )r, = 2r 2 (36a) 



A{ab + cd)(a 2 + c 2 )S + A{a 3 b - c 3 d) cos tp 

+ 2(a 3 d + 3a 2 be + 3ac 2 d + be 3 ) simp 

+ [-4a 3 6 + 6(2a 2 cd + 2abc 2 ) - 4c 3 d]r] = (36b) 



[6a 2 6 2 + 2(a 2 d 2 + Aabcd + b 2 c 2 ) + 6c 2 d 2 ]S 

+ (6a 2 b 2 - 6c 2 d 2 ) cos ip 
+ 6(a 2 bd + ab 2 c + acd 2 + bc 2 d) sin p 

+ 6(-a 2 b 2 + a 2 d 2 + Aabcd + b 2 c 2 - c 2 d 2 )r, = -4r 2 (36c) 



4(a6 + cd)(b 2 + d 2 )S + 4(6 3 a - d 3 c) cos 93+ 

2(6 3 c + 3b 2 ad + 3bd 2 c + ad 3 ) sin p+ 

[-4b 3 a + 6{2d 2 ab + 2cdb 2 ) - M 3 c]r] = (36d) 

(b 2 +d 2 ) 2 S+(b 4 -d 4 )cosp + 2bd(b 2 +d 2 )smp+(-b 4 + 6b 2 d 2 -d 4 )r 1 = 2r 2 (36e) 

along with eq. (29). Notice that the only difference is that the right hand sides 
of eqs.(36a) and (36c) are different from those of eqs. (30a) and (30c). As 
we pointed out in cq. (11), there is another symmetry in parameters: (77, ip — 
tt/2) — > (—r],p = 0). Hence, above analysis is also applicable to the case r/ = 
— 1/4 and tp = 0. In fact the solutions can be reduced to eqs. (13) and (14). 

We plot the kink solutions in Fig. 2 for r\ = 0.15 and p = 0~1.757r. The 
solutions at 77 = 0.35 and p = 0~0.57r and p = 0.5ir~ir are shown in Fig. 3a 
and Fig. 3b respectively. The kinks for 77 < 1/4 (Fig. 2) change smoothly and 
retain their shapes as p moves across n/2. The total change of 9, A8 = 0(s = 
00) — 9(s = —00), is equal to 2ir. This can also be deduced from eq. (28) by 
tracing the variation of 9 with respect to s. On the other hand, the solutions 
for 77 > 1/4 (Fig. 3a and 3b) develop a second kink as p approaches n/2. At 
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ip = 7r/2, the form in eq. (28) is no longer applicable. Eq. (35) has to be 
used and it gives a large kink and a small kinks which are the decedents of the 
connected kinks at ip = it/2 — e where e is an infinitesimal positive number. 
Note also that the total change of 9 of neither the large kink nor the small kink 
is equal to 2tt, but rather the sum of them is. This can also be seen from eq. 
(35). 

We here propose a classical explanation of the situation near ip — tt/2. It 
is related to the minima of the potential. As we have argued in section 2, 
d9 / ds\ s= ± co — and V(0(±oo)) = V m i n . The solutions in eqs. (28) and (35) 
extend from one minimum to another. The kink comes from eq. (28) starts from 
one absolute minimum of the potential, passing through a major peak and ends 
at another absolute minimum. In the special case with the solutions coming from 
eq. (35), there are a large kink and a small kink. Both connect two absolute 
minima but the former passes a major peak of the potential and the latter 
passes a minor peak. The large kink extends from 9 = 2 arctan((6 — a)/(d — c)) 
to 9 = 2 arctan( (b + a)/ (d + c)) and the small kink extend from 6 = 2arctan((6 + 
a)/ (d + c) to 9 = 2n + 2 arctan( (b — a)/{d — c)). It is the emergence of additional 
symmetry in potential which gives rise to two absolute minima in the range of 
2tt which in turn, requires two solutions. We will elaborate more on this in next 
section. 

Finally, we give the form of the bubble solution. Instead of eq. (28), we use 

_ acosh(rs) + b 
H >~ ccosh(r S )+d [6<) 

and follow the same procedure, we are able to obtain the bubble solution, similar 
to that in eq. (22). Its shape is shown in Fig. 4. by the dashed line which 
connects two relative minima of potential. The bubble solution can be found 
only when the relative minima exist. 



0.0.4 4. Discussion and conclusion 

We summarize our results on the ip— r\ phase diagram in Fig. 5. The kink solution 
in eq. (28) can be found in any place except for those vertical lines. On the 
vertical lines, the form in eq. (35) prevails. Note that it corresponds to either 
g(s) = coth(rs) or g(s) = tanh(rs), which in turn, corresponds respectively 
to large kink or small kink. The bubble solution exist in the region above the 
upper dashed line or below the lower dashed line, i.e., the regions where relative 
minima exist. 

The limiting case is also interesting. In the limit r\ — > 0, the solutions 
should be those of an ordinary sine-Gordon equation: 9 = 4 arctan[exp(s)] = 
2 arctan[csc /i(— s)]. It is compatible with the form in eq. (28). In the other 
limit, r/ — > oo, one can write eq. (1) as 
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where 9' = 29. The third term can be treated as a perturbation. The zeroth 
order of the solutions should be those of sine-Gordon equation with sin(2#) : 

9 = 2 arctan[exp(s')] (39) 

Here, s' = 2^/rjs, implying r = 2^/rj in cq. (28) as 77 — > oo. This is shown clearly 
by the Mobius transformation below. 

Consider only the leading order terms of eqs. (30) for large rj 

(a 2 + c 2 ) 2 S + (-a 4 + 6a 2 c 2 - c 4 )rj = 0, (40a) 

A{ab + cd){a 2 + c 2 )S + [-Aa 3 b + 6(2a 2 cd + 2abc 2 ) - Ac 3 d]rj = 0, (40b) 

[6a 2 b 2 +2(a 2 d 2 +4abcd+b 2 c 2 )+6^d 2 }S+6(-a 2 b 2 +a 2 d 2 +Aabcd+b 2 c 2 -c 2 d 2 )rj = 2r 2 

(40c) 

A(ab + cd){b 2 + d 2 )S[-4b 3 a + 6(2d 2 ab + 2cdb 2 ) - M 3 c]rj = (40d) 

(b 2 + d 2 ) 2 S + (-b A + 6b 2 d 2 - d A )r, = 2r 2 . (40e) 

It can be shown easily that none of a, b, c and d can be 0. Hence we can set 
a = pc and b — qd where p and q are just two ratio parameters. With the 
condition eq. (29), it can be shown further that q = —p = ±1 and \cd\ = l/\/2. 
Now let S = —i] + 5 where 5 = 0(rp) and substitute it into eq. (40e), one find 
that S = 2r 2 /d i . This implies that b = ±d = O(tfrj) and a = ±c = 0(1/ ^fj). 
We thus have shown that for finite rj, solutions of the form of eq. (28) always 
exist (except for (rj < 1/4, ip = 0) and (77 > 1/4, tp = n/2)) and their parameters 
a, b, c and d vary smoothly with rj. Hence, the solution for finite rj has the same 
form as that for rj = 0. 

The only places where there are phase transition are the ends of the vertical 
lines in Fig. 5, i.e., the points (ip = nir, rj = —1/4) and (ip = (n+l/2)n, rj = 1/4). 
Here, indeed the form of the solutions changes from that in eq. (28) into that in 
eq. (35) when the absolute value of rj increases. Classically, this is a second-order 
phase transition. Its quantum fluctuation has also been well-studied [12- 15]. 

There is another aspect we would like to investigate and that is varying ip 
across ip — n/2 for a fixed rj. We plot in Figure 6a and 6b the energy as a function 
of ip with rj = 0.15 and 77 = 0.35 respectively. One can immediately notice the 
behavior of energy near ip = tt/2 + nir. For rj = 0.15, the slope change is large 
but still smooth. For rj = 0.35, though the energy remains continuous, the slope 
does not. This indicates that when \rj\ < 1/4 there is smooth crossover. But 
when 1 77 1 > 1/4 there is a second-order phase transition from the energy point 
of view. The solutions also show different behavior. In Fig. 2, the kinks vary 
smoothly across the point ip = tt/2 at 77 = 0.15. For the solutions in Figs. 3 with 
77 = 0.35, one finds that near ip rts n/2, the shapes of the solutions are different 
from those at ip = n/2. The solutions are combination of two kinks though their 
forms are still that of eq. (28), i.e. , g(x) = sinh(rx). As a result the range of 
variation of 9 is still 2n. and the topological charge is unity. When <p = n/2 the 
form of the kinks is that of eq. (35), i.e., g(x) = tanh(rx) or g(x) = coth(rs), 
the large or small kink solutions. It is the sum of ranges of variation of 9 of the 
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two kinks which is equal to 2ir, and thus the deconfinement. It is due to the 
emergence of the symmetry that both 9 and ir — 9 are solutions. On the other 
hand, the solutions at <p = ir/2±e where e is a infinitesimal number, are very 
similar. Our numerical results also show that the coefficients a, b, c, d and S 
vary smoothly across the point ip — n/2, (not including the point ip = tt/2.) 
Thus, the point ip = n/2 for r\ > 1/4 is actually a singular point. 

In this work, we used the method of " Mobius transformation" to solve the 
asymmetric DSGE. This method transformed the DSGE into a set of algebraic 
equations. Thus we are able to find the forms of the solutions for all the region 
on the ip — r\ plane. The resulting forms of our solutions can serve as the basis of 
various methods, such as form factor perturbation theory, semi-classical method 
or a truncated conformal space approach to study quantum fluctuation. 

This work is supported in part by NSC of Taiwan, ROC under the contract 
number NSC 95-21 12-M-002-048-MY3. One of the authors, C. D. Hu would like 
to thank Chern Chyh-Hong for inspiring discussion. 

0.0.5 Appendix: Introduction to Mobius transformation 

By Jacobi elliptic function theory, one can transform the following equation: 

if/ = p(s) = A(f fa)(f fa)(f f 2 )(f fa), (A-l) 

where ip(s) is a polynomial of s to the three or four power and fa, fa, fa and 
fa are the roots, into the standard form, i.e., only terms with even powers are 
present. The "Mobius transformation" has the form 



and so does every root 



/, - (A-3) 
cQ + d 



If we take special values of a, 6, c and d, we can obtain the form 

r = fai.fl - /o) * C - faifa - fa) /a a\ 

1 (h fo)C (fa fa) 1 " ' 

and eq. (A-l) becomes 

( S )2 = B{ ' ~ P2){ " ~ M - (A " 5) 

We set 

> _ fa - fa fa - fa _ Pi -00 fo- fh (A ) 

fa- fa fa- fa 01-02 03-00' 1 " ' 



so that eq. (A-5) becomes 



(f/ = B^-l)(X,-l) (A-7) 
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here B = A(fs — /0X/2 — fi)- This is the standard form. But we can do further 
transform by setting <; = £ 2 . So cq. (A-7) becomes the differential equation of 
Jacobi elliptic function: 

(f/ = j(e-i)(M 2 -i)- (a-8) 

Instead eqs. (A. 4) and (A. 8), we can take different transformation by map- 
ping roots into 1, — 1, 1/fc, and —1/k. So 

x _ h - fo h - h _ I - 1 it + 1 _ , 1 - ^ 2 /a Q N 

/1-/2/3-/0 ~ l + l^-l'h + k' (A_9) 



and eq. (A-l) becomes 



(g) 2 = £'(£ 2 -l)(fc 2 e 2 -l) (A-10) 



There are more details to be solved such as how to deal with degenerate roots 
and how to further transform A so that it is real and in the range (0, 1) in order 
that eq. (A-8) and (A-10) are compatible with the Jacobi elliptic differential 
equation. But these are beyond the scope of this work. 
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Figure cations: Fig. 1 Potential energy v.s. 9 for ip = and (a) 77 = —0.35, 
(b) r] = 0.15 and (c) 77 = 0.35. 

Fig. 2 Solutions form the form of eq. (28) for 77 = 0.15 and p = 0~7n/8. 

Fig. 3a Solutions from the form of eq. (28) for 77 = 0.35 and <p> — 0~7r/2. At 
p = it/2, the solutions have the form of eq. (35) which gives large and small 
kinks. 

Fig. 3b Solutions from the form of eq. (28) for 77 = 0.35 and ip = ir/2~ir. At 
p = 7r/2, the solutions have the form of eq. (35) which gives large and small 
kinks. 

Fig 4. kink solution (in dashed line) and bubble solution (in solid line) for 
77 = 0.35 and p = 0. 9 is in unit of ir. 

Fig. 5 Phase diagram on p — 77 plane. See text for explanation. 
Fig. 6 The energy of the kink versus <p, (a) 77 = 0.15, (b) 77 = 0.35. 
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Tabic 1 List of a, b, c, d, r, S for eq. (28) ionp = 





a 


b 


c 


d 


r 


S 


-0.25 


X 


X 


X 


X 


X 


0.75 


-0.15 





0.795271 


1.25743 





±0.532456 


0.85 








1 


1 





±1 


1 


0.15 





1.12468 


0.88914 





±1.26491 


1.15 


0.25 





1.18921 


0.840896 





±1.41421 


1.25 


0.35 





1.24467 


0.803428 





±1.54919 


1.35 


0.45 





1.29357 


0.773055 





1.67332 


1.45 



Table 2 List of a, b, c, d, r, S for eq. (28) ionp = ir/4 



v\ 


a 


b 


c 


d 


r 


5" 


-0.25 


-0.466708 


0.836099 


0.712245 


0.866693 


±1.26959 


1.10092 


-0.15 


-0.460968 


0.828478 


0.7977 


0.735676 


±1.1254 


1.04093 





-0.382683 


0.92388 


0.92388 


0.382683 


±1 


1 


0.15 


-0.238106 


1.10602 


0.890012 


0.0656207 


±1.1254 


1.04093 


0.25 


-0.173621 


1.20406 


0.833646 


-0.0216328 


±1.26959 


1.10092 


0.35 


-0.13299 


1.2828 


0.785775 


-0.0600972 


±1.41323 


1.17461 


0.45 


-0.106325 


1.34854 


0.74763 


-0.0772086 


±1.54803 


1.25625 



Table 3 List of a, b, c, d, r, S for eq. (28) ionp = n/2 



v\ 


a 


b 


c 


d 


r 


S 


-0.25 


-0.594604 


0.840896 


0.594604 


0.840896 


±1.41421 


1.25 


-0.15 


-0.628717 


0.795271 


0.628717 


0.795271 


±1.26491 


1.15 





-0.707107 


0.707107 


0.707107 


0.707107 


±1 


1 


0.15 


-0.88914 


0.562341 


0.88914 


0.562341 


±0.632456 


0.85 


0.25 


X 


X 


X 


X 


X 


X 


0.35 


X 


X 


X 


X 


X 


X 


0.45 


X 


X 


X 


X 


X 


X 
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1.5 




-1.5 1 1 1 1 1 

-2-1012 

6 (in unit of n) 



17 




6 (in unit of n) 



18 



1.5 




6 (in unit of n) 
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20 



8 



r|=0.35 




21 




22 




23 



1.00 
0.75 f- 
0.50 [■ 
0.25 
0.00 
-0.25 k 



-1.00 













/'"""~\ 




/ 


s / \ ' 

x^ ^ ' x^ ^ ^ 






1 


1 


1 


i 



-1.00 -0.75 -0.50 -0.25 0.00 0.25 0.50 0.75 1.00 

q>(7c) 



24 




25 




5 - 1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 18 

Adiabatical phase( 71/8) 
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